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1. Introduction 

Models of the Nambu-Jona-Lasinio (NJL) type have a long history and have 
been extensively used to describe the dynamics of lightest hadr qns and t he ther- 
modynamic properties of excited matter (see the review articles IMISEEI). i n the 
"classical" versions such schematic models incorporate the chiral symmetry of two- 
flavor QCD and its spontaneous breakdown at temperatures below the critical one, 
T < T c . They offer a simple and practical illustration of the basic mechanisms that 
drive the spontaneous breaking of chiral symmetry, a key feature of QCD in its 
low-temperature, low-density phase. However, in spite of their widespread use, the 
NJL models suffer from the major shortcoming that the reduction to global color 
symmetry prevents quark confinement. 

In the Polyakov-loop-extended NJL (PNJL) model ISHEl, the quarks are cou- 
pled simultaneously to the chiral condensate, to be an order parameter of the chi- 
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ral symmetry breaking, and to a homogeneous gauge field representing Polyakov 
loop dynamics, which serves as an order parameter for the transition from the low- 
temperature, symmetric, confined phase to the high-temperature, deconfined phase. 
Th e m odel has proven successful in reproducing lattice data on QCD thermodynam- 

icsHni 

In this paper we confront general properties of 7r and cr mesons as well as ther- 
modynamics at finite temperature T and baryon chemical potential (i calculated 
within the two-flavor NJL model with those of the PNJL one. 



1.1. Nambu-Jona-Lasinio model 

To describe the coupling between quarks and the chiral condensate in the scalar- 
pseudoscalar sectors, the two- flavor NJL model H] 6 | 13 | 14 | ^ g usec j w jth the following 
Lagrangian density 

£njl =q(ip- mo) q + G (qq) 2 + (qijsrq) , (1) 

where G is the coupling constant, f is the vector of Pauli matrices in flavor space, q 
and q are the quark fields (color and flavor indices are suppressed), too is the diagonal 
matrix of the current quarks masses, mo = diag (to", m"), to" = to" = mo. 

The grand canonical thermodynamic potential can be obtained from this La- 
grangian in systematic approximations. In the mean-field approximation it has the 
form 



^njl — G(qq) 2 + Cl q , (2) 

with 

n q = -2N c N f J -0^E P - 2N c N f T J [in N + (E P ) + lnJ\T(£p)] , (3) 

where N+(E p ) = l + er^-^ &nAN-{E p ) = l + e '^ E "+^ with E p = ^/p 2 + to 2 
and the inverse temperature (3 — 1/T. 



1.2. Nambu-Jona-Lasinio model with Polyakov-loop 

The deconfinement in pure SU(N C ) gauge theory can be simulated by intro ducing an 
effective potential for a complex Polyakov loop field. The PNJL LagrangianDSIIHHlH 
is 

C PN j h = q(i^-mo)q + G[(qq) 2 + (qi^Tq) 2 ] -U ($[A],*[A];T) . (4) 

Here, the notation is the same as in Eq. ([T]). 

The quark fields are coupled to the gauge field A^ through the covariant deriva- 
tive = - iA*. The gauge field is A 1 = 6gA° = -i5%A 4 (the Polyakov gauge). 
The field 3> is determined by the trace of the Polyakov loop L(x) 

*[A] = -^Tr c L(x) , 
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where L(x) — V exp 



The gauge sector of the Lagrangian density @ is described by an effective po- 
tential U ($[A], T) fitted to the lattice QCD simulation results in pure SU(3) 
gauge theory at finite T l l^ l -^ l with 
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The parameters of the effective potential and © are summarized in Table [T] 



Table 1. Parameters of the effective po- 
tential U[A\. 



do 


«i 


a.2 


03 h 


bi 


6.75 


-1.95 


2.625 


-7.44 0.75 


7.5 



The parameter To in general depends on the number of active flavors and the 
chemical potental^^. In the present work we use To = 270 MeV as in^l 

The thermodynamic potential for the PNJL model in the mean-field approxi- 
mation is given by the following equation ^\ 



n($, m, T,fj,)=U ($ , !>; T) + G(qq) 2 + Sl q , (7) 
where (in analogy with {2])) 

n q = -2N c N f J ^E p - 2N f T J j^L [lnN+(E p ) + \nN^(E p )] . (8) 
Here, E p — \/p 2 + m 2 is the quark energy, E^ = E p ^f fi, and 



N£(E P ) = 


1 - 


-3 




-^ E t)e-^ E t +e-^ E t' 


(9) 




1 - 


1-3 


%- 




(10) 



Since NJL-type models are nonrenormalizable it is necessary to intro duc e a 
regularization, e.g., by a cutoff A in the momentum integration. Following we 
use in this study the three-dimensional momentum cutoff for vacuum terms and 
extend this integration till infinity for finite temperatures. A comprehensive study 
of the differences between the two regularization procedures (with and wit hou t 
cutoff on the quark momentum states at finite temperature) was performed in^l 
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2. Quarks and light mesons in NJL and PNJL models 

In the mean-field approximation, we can obtain the constituent quark mass to from 
the condition that the thermodynamical potential (Eqs. © and ([7]), resp.) shall have 
a minimum with respect to va rying t his parameter, dQ/dm = 0. This condition is 
equivalent to the gap equation I I ' > 



m = toq — 2G (qq) 



(11) 



where the quark condensate is defined as (qq) — dfl/dmo. For the mass gap equation 
of both models we get 



to = m + 8GN c N f / -^slT [l - /+ - /" 



with 



(12) 



f + = (l + e^)-\ 
f- = (l + e^)- 1 



(13) 
(14) 



for the NJL model, and 



$ + 2$e 



-PE- e -f)E- 



-3/3B+ 
-3/3B- 



(15) 
(16) 



for the PNJL model. Moreover, for PNJL calculations we should hnd the values of 
<I> and $ by minimizing £1 with respect $ and $ at given T and /i. One should 
note, that if $ — > 1, the expressions Eqs. (fl"5|) . (fT6 |) reduce to the standard NJL 
model Eqs. (JT3]) and (fT4"|). 

For a self-consistent description of the particle spectrum in the mean-field ap- 
proximation, the meson correlations have to be taken into consideration. These cor- 
relations are related to the polarization operator of constituent fields. For scalar 
and pse udoscala r particles the polarization operators are represented by loop- 
integralsHHSm. 



nf b P (P 2 ) = / [i^r a S{p + P^ 5 r b S(p)] , (17) 

nff (P 2 ) = J ^jTr [S(p + P)S(p)} , (18) 

where the operation Tr is taken over Dirac, flavor and color indices of quark fields. 

From point of view of the polarization operators, The pseudoscalar (ir) and 
scalar (a) meson masses can be defined by the condition that for P 2 — M 2 (M 2 ) 
the corresponding polarization operator II PP (M^) (H SS (M 2 )), leads to a bound 
state pole in the corresponding meson correlation function l^A For mesons at rest 
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(P = 0) in the medium, these conditions correspond to the equations 

1 + 8G WlW^? (w+ - r> = ' (19) 

i + «e^/|^|^(W + -r} = o. (20) 

In order to solve Eqs. (fTTj) , (fT9|) and (|20|) , a set of model parameters has to 
be determined: the cutoff parameter A, the current quark mass m (in chiral limit 
mo = 0) and the coupling constant G. These parameters are fixed at T = to 
reproduce physical quantities: the pion mass — 0.139 GeV, the pion decay 
constant = 0.092 GeV and the quark condensate (qq) 1/3 = -250 MeV. The 
obtained parameters are shown in Table [2j 



Table 2. The set of model parameters reproducing ob- 
servable quantities (in brackets) and the chiral condensate 
(qq) 1 / 3 = -250 MeV. 



m [MeV] 


A [GeV] 


G [GeV]" 2 


F n [GeV] Mir [GeV] 


5.5 


0.639 


5.227 


(0.092) (0.139) 



Solutions of the gap-equation (|11[) and Eqs. (fl"9|) . (|20|) at nonzero T are presented 
in Fig. [TJ The temperature is normalized to the Mott temperature, which is defined 
from the condition M^(TMott) = 2m g (TMott)- In the PNJL model Tm u — 0.27 GeV 
and in the NJL model Tm u — 0.208 GeV for our parameters. The modification of 

0.7 
0.6 
0.5 
> 0.4 

0) 

CD 

- 0.3 
0.2 
0.1 
0.0 

0.0 0.5 1.0 

T/T M„„ 

Fig. 1. Temperature dependence of the masses m q , and M a at /i = GeV. Results for PNJL 
and NJL models are given by the solid and dashed lines, respectively. 

the quasiparticle properties is clearly seen in this figure. Up to the Mott temper- 
ature T]yiott, the a mass practically follows the behaviour of 2m q (T) with a drop 
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towards the pion mass signalling chiral symmetry restoration. At T/Tmou > 1 the 
masses of chiral partners become equal to each other, M a » M n , and then both 
masses increase with temperature. Below the Mott temperature, the pion mass re- 
mains practically constant. The transition region from the phase with broken chiral 
symmetry (m q (T) ~ m q (0)) to chirally symetric phase (m q (T) ~ 0) is much nar- 
rower in the PNJL case when compared to the NJL model. For a recent discussion 
of this issue within the nonlocal PNJL model, see ^1 

3. Thermodynamics of NJL and PNJL models 

The thermodynamics of particles is described in terms of the grand canonical en- 
semble which is related with the Hamiltonian H as follows: 



e 



-0VQ _ Tr e -fHH-^N)^ (21) 



where N is the particle number operator and the operator Tr is taken over momenta 
as well as color, flavor and Dirac indices. If is known, the basic thermodynamic 
quantities - the pressure, the energy and entropy densities, the density of quarks 
number and heat conductivity - can be defined as follows: 

(22) 



(23) 



T ( ds 



'cKT 

-p + Ts + ^n, (24) 

(25) 

v yar 1 ' (26) 

The thermodynamic potential in equilibrium corresponds to a global minimum with 
respect to variations of the order parameter(s) 

dn(T,n,m) _ Q; d 2 n(T,^m) ^ Q 

dm dm 2 
All these relations describe thermodynamics of the system. For the considered 
models the thermodynamic potentials are defined from Eqs. ^ and ([7]). From these 
equations we can read off the vacuum part 

a m ^^p^-2N c N,J^E,. (28) 

This quantity does not vanish at T — )• and fi — > 0. Therefore, in order to obtain the 
physical thermodynamical potential which corresponds to vanishing pressure and 
energy density at (T, /j) = (0, 0), one has to renormalize the thermodynamic poten- 
tial by subtracting its vacuum expression (|28j) . This corresponds to the following 
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defintion of the physical pressure 

p p(T, fi, m) — p(0, 0,m) 



j>4 - T 4 ■ (29) 




Fig. 2. The temperature dependence of the reduced pressure and energy density within the PNJL 
model for fi = in two schemes of regularization A =0.639 and A — > oo. Dotted lines are appropri- 
ate results for the NJL model. Lattice data points for Nf = 2 at fj, =0 are from Ref. [27]. Circles, 
squares and diamonds correspond to calculations at Nt =6 with the mass ratio of the pseudoscalar 
to vector meson mps/my =0.65, 0.70 and 0.75, respectively. 



Within the PNJL model with A — > oo the reduced pressure and energy density 
exhibit reasonable behavior consistent with the recent lattice QCD results for the 
vanishing chemical potential (see Fig. ^ keeping in mind that the mps/mv ratio 
in lattice calculations is still far from that for physical masse s rnp s I m v ~ 0-2- One 
should note that both models have the cutoff parameter A 03. But the in tegrals 
containing the logarithm in Eq. {2} and Eq. ([7]) are both convergent 1 ^ 1 25 1 j n our 
work we calculated these integrals with A — > oo. It leads to the flattening of pressure 
at high temperature. However, most integrals for the PNJL model are convergent 
too. This was the reason to consider the thermodynamic functions for A — > oo. 
It was supposed that with increasing temperature the pressure has to reach the 
Stefan-Boltzmann limit ^2H, which in the PNJL model is defined as 

ji = ( N c ~ 1)45 + NcN f — * 4.053. (30) 

If the regularization A = 0.639 is used, the T-behaviour of the thermodynamic 
quantities considered is roughly the same while their absolute values are noticeably 
lower, being far from the Stefan-Boltzmann limit. In the NJL model both, p/T 4 
and e/T 4 are not only underestimated due to the missing gluon contribution, but 
also essentially shifted toward lower temperatures because of the lack of a confining 
mechanism for the dynamical quark degrees of freedom. 

Within NJL-like models there are several characteristic temperatures. The pa- 
rameter To entering the effective potential (j6]) of the PNJL model has been noted 
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above. Three other scales are the pseudo-critical temperature for chiral crossover 
T x defined by the maximum of d{qq)/DT, the pseudo-critical temperature for the 
deconfinement crossover T p which can be found from the maximum of 9$/9T, and 
T c define d for PNJL model as the average of two transition temperatures T x and 
j^| 10 | 15 ] rpjjg temperature dependence of the order parameters for the chiral (((?(?)) 
and deconfinement ($) phase transitions, are shown in Fig. ([3]). The chiral con- 
densate decreases and the Polyakov loop potential increases with T, demonstrating 
closeness of the pseudo-critical temperatures T x and T p at fi = (see also Ta- 
ble [3]). For 7r-mesons, the Mott temperature T Mott is provided by the condition 



1.25 






1.00 


<qq>/<qq> T _ 








0.75 






0.50 






0.25 
0.00 







0.2 0.3 
T [GeV] 



Fig. 3. Temperature dependence of the chiral condensate and Polyakov loop potential at fi = 
GeV within the PNJL model 

(TMott) = 2m 9 (TMott) and similarly the a meso n dissociation temperature T% is 
given by the equation M<j(T%) = 2M„(Tg) I 19 ! 21 ! All these quantities obtained at 
/i = are presented in Table [3] 



Table 3. Characteristic temperatures in NJL and PNJL 
models. 





To 


T x 


T P 


T c 


^Mott 


1 d 


NJL 




0.192 




0.192 


0.207 


0.185 


PNJL 


0.27 


0.249 


0.258 


0.253 


0.27 


0.259 



Extending our study of the (pseudo-)critical temperatures to nonzero chemical 
potential ^i, we obtain phase diagrams in the T — fx plane shown in Fig.|4j The chiral 
transition line, determined by^ 

"_L_ dCl a 
AG + 



dm 2 







m— 



in both NJL and PNJL models, is a monotonously decreasing function of the chem- 
ical potential. In the limiting chirally symmetric case corresponding to m = the 
T x at large fi is higher than those for finite mass but both temperatures coincide 
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when n — >Q for both the NJL and PNJL models. For the case /i ^0 both mod- 
els show the critical end point at the temperature Tcep below which the chiral 
phase transition is of first order. At this poi nt (Tcep, Hcep) the phase transition 
changes from first order to crossover HHHHM] p or ^ e chjj-gjiy symmetric case in the 
PNJL model the first order phase transition ends at a tricritical point above which, 
for T > Ttcp, the chiral transition is of the second order. In the NJL model 
the topology of the phase diagram is the same as in the PNJL case, only Tcep 
and Ttcp are situated at higher temperatures. In accordance with other calcula- 
tions, the temperature of the tricritical point is above that of the critical end point. 
Within the PNJL model the positions of critical points are (Tcep, I^cp) =(95,320) 



NJL 



0.25 
0.20 
0.15 
0.10 
0.05 



0.00 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ ^- 1 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



n[GeV] 



CD 



0.30 






PNJL 


0.25 
0.20 








0.15 




TCPN 










\ \ CEP ■ 


0.10 








0.05 








0.00 









0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 
H [GeV] 



Fig. 4. Phase diagrams of NJL (left panel) and PNJL (right panel) models. Solid lines correspond 
to the first order phase transition, dashed lines are crossover and dotted lines are the boundary of 
the second order phase transition. 



and (Ttcp, Htcp) =(160,265) MeV. These numbers are quite close to those in^^l 
for the set with similar parameter values (set B, the case I). One should empha- 
size that critical properties of observables are significantly influenced by the chosen 
parameter set and regularization procedure as was demonstrated in 12-21 



4. Conclusion 

We have compared the thermodynamics of NJL and PNJL models. In agreement 
with previous results, it is shown that the inclusion of coupling between chiral 
symmetry and deconfinement essentially improves the description of thermodynamic 
bulk properties of the medium. The models qualitatively reproduce both, 7r and a 
meson properties in hot, dense quark matter and the rich and complicated phase 
structure of this medium. Effects of the Polyakov loop move the CEP to higher T 
and lower fi than in the NJL case ^^1. The position of the calculated CEP in the 
T — fx plane is still far from the predictions of lattice QCD and empirical analysis. 
The further elaboration of the presented models may include color superconducting 
phases and nonlocality of the interaction - as well as effects beyond the meanfield 
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